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We investigate the nucleosynthesis constraint on Lorentz invariance violation in the neutrino sector 
which influences the formation of light elements by altering the energy density of the Universe and 
weak reaction rates prior to and during the big-bang nucleosynthesis epoch. We derive the weak 
reaction rates in the Lorentz-violating extension of the standard model. Using measurements of 
the primordial helium-4 and deuterium abundances, we give a tighter constraint on the deformed 
parameter than that derived from measurements of the cosmic microwave background anisotropies. 
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Neutrino oscillation experiments have shown that there 
are small but non-zero mass squared differences between 
three neutrino mass eigenstates (see Ref. [1] and refer- 
ence therein). However, neutrino oscillations cannot pro- 
vide absolute masses for neutrinos. Cosmology provides 
a promising way to constrain the total mass of neutrinos 
by the gravitational effect of massive neutrinos on the ex- 
pansion history near the epoch of matter-radiation equal- 
ity [2] and on the formation of large-scale structures in 
the Universe [3] (see also Ref. [4] for a review). Recently, 
a 30 detection of non-zero neutrino masses is reported 
using new measurements of the cosmic microwave back- 
ground (CMB) anisotropies from the south pole telescope 
and Wilkinson microwave anisotropy probe (WMAP), 
in combination with low-redshift measurements of the 
Hubble constant, baryon acoustic oscillation feature and 
Sunyaev-Zel’dovich selected galaxy clusters [5]. 


These observations establish the existence of physics 
beyond the standard model of particle physics. Another 
possible signal of new physics is violation of Lorentz sym- 
metry. The possibilities of Lorentz invariance violation 
were considered in string theory [6], standard model ex- 
tension [7], quantum gravity [8], loop gravity [9], non- 
commutative field theory [10], and doubly special relativ- 
ity theory [11]. Searches for Lorentz invariant violation 
with neutrinos have been performed with a wide range 
of systems [12]. Although present experiments confirm 
Lorentz invariance to a good precision, it can be broken 
in the early Universe when energies approach the Planck 
scale. Cosmological observations provide a possibility to 
test such a symmetry at high energies. 


Recently, measurements of the CMB power spectrum 
were used to probe Lorentz invariant violation in the neu- 
trino sector [13]. Lorentz invariant violation affects not 
only the evolution of the cosmological background but 
also the behavior of the neutrino perturbations. The for- 
mer alters the heights of the first and second peaks in the 
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CMB power spectrum, while the latter modifies the shape 
of the CMB power spectrum. These two effects can be 
distinguished from a change in the total mass of neutrinos 
or in the effective number of neutrinos. The seven-year 
WMAP data in combination with lower-redshif measure- 
ments of the expansion rate were used to put constraints 
on the Lorentz-violating term. However, the resulting 
constraints suffer from a strong correlation between the 
Lorentz-violating term and the dark matter density pa- 
rameter [13]. 

In this letter, we use current big-bang nucleosynthe- 
sis (BBN) data to constrain Lorentz invariance violence 
in the neutrino sector. There are two effects of Lorentz 
invariant violation on BBN. The first is a correction to 
the weak reaction rate in the Lorentz-violating standard 
model extension, which governs the neutron-to-proton 
ratio at the onset of BBN. The second is a change in 
the total energy density of the Universe. Since the abun- 
dances of the light elements produced during BBN de- 
pend on the competition between the expansion rate of 
the Universe and the nuclear and weak reaction rates, the 
BBN predictions depend on the Lorentz-violating term. 
In particular the BBN-predicted abundance of helium-4 
is very sensitive to the deformation parameter. 

We focus on Lorentz invariance violence only in the 
neutrino sector and consider the following deformed dis- 
persion relation 


E? =m? +p? + €p’, (1) 


where F is the neutrino energy, m the neutrino mass, 
p = (ppi)! the magnitude of the 3-momentum, and 
€ the deformation parameter characterizing the size of 
Lorentz invariance violation. The dispersion relation im- 
plies that there are departures from Lorentz invariance 
in the neutrino sector if € #0. Such a deformed disper- 
sion relation was constructed in the framework of con- 
ventional quantum field theory [14] and derived in the 
Lorentz-violating extension of the standard model [15]. 
Here we have to point out that the dispersion rela- 
tion for neutrinos given in (1) is not very general. It ne- 
glects neutrino oscillations, possible species dependence, 
anisotropies associated with violations of rotation sym- 
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metry, and CPT violation. As shown recently by Kost- 
elecky and Mewes, all of these are possible [12]. The 
model considered in this paper is one of many possible 
Lorentz-violating theories. 

The number density n, and energy density p, for mas- 
sive neutrinos with (1) are given by [13] 


3 
Ny = s | ast), (2) 


dp 
v = gv HE 
a J Cr) 
where gy = 2 is the number of spin degrees of freedom. 


The phase space distribution for neutrinos is the Fermi- 
Dirac distribution 


f(E) = [1 + exp(E/T,)|*, (4) 
where T, is the neutrino temperature. Thus the e 
and energy density = be written as n, = (1+€)~3/2n (9) 


and py = (1+€)7 3/2) where n® and p® are the stan- 
dard number and eee density, respectively. Increasing 
€ decreases both the number and energy density. The 
former leads to a reduced rate of the weak reaction prior 
to and during the BBN epoch since the reaction rate is 
proportional to the neutrino number density, while the 
latter results in a reduced expansion rate of the Uni- 
verse. Therefore, Lorentz invariant violation affects the 
nucleosynthesis of light elements. 

We turn now our attention to the details of the compu- 
tation of the weak reaction rate with Lorentz invariance 
violation in the neutrino sector. At early times when the 
temperature of the Universe was T ~ 100 MeV, the num- 
ber and energy density were dominated by relativistic 
particles: electrons, positrons, neutrinos, antineutrinos 
and photons. All of the particles were kept in thermal 
equilibrium by the weak reactions 


Ve+tn & pte, (5) 
erin & p+, (6) 
no ptit+e. (7) 


When the expansion rate of the Universe exceeds the 
reaction rate for n 4+ p processes, the baryons become 
uncoupled from the leptons. At the time the neutron- 
to-proton ratio is frozen, which largely determines the 
primordial helium mass fraction. To estimate the neu- 
tron abundance at the onset of BBN one has to compute 
the reaction rate. As an example, let us consider the pro- 
cess Ve +n > p+e-. The differential reaction rate per 
incident nucleon is 


_ MP dp, dpe 
du = Dy Samp a” arp, T Fe) 


Mp — Ee), (8) 


spins 


276 (Ey + Mn — 


where |M]? is the squared matrix element, to be summed 
over initial and final state spins, Mn and mp the neutron 


and proton mass, respectively, (Ee, pe) the electron four- 
momentum, and fe denotes the Fermi-Dirac statistical 
distribution for electron. The process (5) involves the 
gauge boson W as mediator. At tree level, one has 


M= - Up Yp(Cv — Cays)Untie(y* + ch” 
where Gp is the Fermi coupling constant, c,, are the 
coefficients for Lorentz violation, and Cy, Ca are the 
vector and axial coupling of the nucleon. Here the coeffi- 
cients Cuy are defined to be traceless and isotropic. After 
integrating (8) the reaction rate is 


3 (Cy — Ca) 


A (0) 
4(GE +302) ae 


1-s€- gated 


where w( is the standard reaction rate per incident nu- 
cleon derived in [16]. The first factor on the right-hand 
side of (10) arises from the neutrino propagator and the 
evW coupling in the Lorentz-violating extension of the 
standard model [7] and the second factor from the sta- 
tistical distribution for neutrinos (more general Lorentz- 
violating corrections involving electrons, neutrinos, neu- 
trons and protons were discussed in [17]). At tree level, 
the differential reaction rates for the other five processes 
in (5)-(7) can be simply derived from (8) by properly 
changing the statistical factors and the delta function 
determined by the energy conservation for each reaction. 


FIG. 1: The “He mass fraction (solid curve) and D/H abun- 
dance ratio (dashed curve) from BBN theory as a function 
of the deformation parameter € for Qh? = 0.0213 (the up- 
per panel) and 0.0224 (the lower panel). The vertical line 
corresponds to Lorentz invariance. 


Ww) (1 — 75) uv (9) 
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Therefore, the corrections to the conversion rate of neu- 
tron into proton and its inverse are the same as in (10). 

From (10) we see that increasing € reduces the reac- 
tion rate, and thus the weak reactions freeze out at earlier 
time, corresponding to a higher freeze-out temperature. 
This leads to a larger neutron-to-baryon ratio at the onset 
of BBN and thus a larger abundance of primordial “He 
production. On the other hand, increasing € also reduces 
the expansion rate of the Universe due to a decrease of 
the energy density, which means the weak reactions freeze 
out at later time without corrections to the reaction rate 
induced by the deformed parameter. This therefore re- 
sults in a lower helium-4 abundance. These two effects 
play opposite roles in the BBN prediction for the helium- 
4 abundance. The abundances of the other light nuclides 
weakly depend on € by changing the neutron-to-proton 
ratio and the expansion rate. 

Considering these corrections to both the reaction rate 
and the expansion rate, we now estimate the freeze- 
out temperature, Tj, determined by equating the ex- 
pansion rate with weak reaction rate. In the Friedman- 
Robertson-Walker Universe, the expansion rate obeys 
H? = 8nGp/3 where p x T* at early times. Thus, we 
have H œ (1—0.75£)T?. Since the standard reaction rate 
in Eq. (10) is roughly given by w œ T’ [18], we have 
w x (1—1.80€)T°. Setting H ~ 4w since the free-neutron 
decay process and its inverse are neglected at the BBN 
epoch, one derives the freeze-out temperature 


Ty ~ (1+0.35€)T}”, (11) 


where T°) is the standard one. For a large €, the weak 
reactions freeze out at a higher temperature. This im- 
plies that effects caused by changing the reaction rate 
dominate over those by changing the expansion rate due 
to the Lorentz invariance violation in the neutrino sector. 

In order to calculate the abundances of light elements 
produced during BBN, we modified the publicly available 
PArthENoPE code [19] to appropriately incorporate the 
Lorentz-violating term in the neutrino sector. Figure 1 
shows the “He mass fraction and D/H abundance as a 
function of € for Qh? = 0.0213 (the upper panel) and 
0.0224 (the lower panel). Both Y, and D/H increase as 
€ increases since the effect of changing the reaction rate 
play a leading role. Moreover, the dependence of Y, on 
€ is much larger, relative to its observational uncertain- 
ties, than that of D/H. Therefore, the primordial helium- 
4 abundance can provide a sensitive probe of neutrino 
physics with Lorentz invariance violation. 

Assuming that there are three types of neutrinos with 
vanishing chemical potentials in the Universe, the BBN- 
predicted primordial abundances depend on only two pa- 
rameters: 0,h? and £. As shown in Figure 1, the abun- 
dance of deuterium is more sensitive to the baryon energy 
density parameter but less sensitive to the deformation 
parameter while that of helium-4 is more sensitive to € 
but less sensitive to Q,h?. We use the observed primor- 
dial abundances of “He and D in combination to constrain 
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FIG. 2: Two-dimensional joint marginalized constraints (68% 
and 95% confidence level) on the deformation parameter € and 
physical baryon density Neh? from measurements of Yp and 
D/H. The dashed line corresponds to Lorentz invariance. 


these two parameters based on the likelihood function 


(Yp — 0.2565)? 
0.0062 


(log[D/H] + 4.55)? 


ine = 
ne 0.032 


(12) 
Here we adopt the estimate of the primordial helium 
mass fraction, Yp = 0.2565 + 0.0060, derived in [20] using 
Monte Carlo methods to solve simultaneously for many 
possible systematic effects, based on 93 spectra of 86 low- 
metallicity extragalactic HII regions. While some have 
employed a selected subset of these data for more de- 
tailed analyses, the sources and magnitudes of system- 
atic errors have rarely been addressed. The measure- 
ment uncertainty in Y, is currently dominated by sys- 
tematic errors. For the primordial deuterium abundance, 
we use the value of log[D/H] = —4.55 + 0.03 obtained 
in [21] from measurements of the absorption lines of seven 
quasars at high redshifts in low-metallicity hydrogen-rich 
clouds with low internal velocity dispersions. Besides 
4He and D, *He and ‘Li are the other two nuclides pre- 
dicted in measurable quantities by BBN. Since their post- 
BBN evolutions are complicated and their measurements 
suffer from systematical uncertainties which are difficult 
to quantify (for helium-3) or are poorly understood (for 
lithium), the observed ?He and ‘Li do not provide a reli- 
able probe of BBN at present, as discussed in [22]. Thus, 
we do not include them in our constraints. 

The He abundance is used to provide a constrain on 
the deformation parameter while the D abundance is used 
to provide a constrain on the baryon density parameter. 
Using the combination of the “He and D data, we find 
E€ = 0.036 + 0.023 and Q,h? = 0.0213 + 0.0009 (68% 
confidence level). This estimated value of the deforma- 
tion parameter is consistent with Lorentz invariant € = 0 
within 95% confidence level. Compared to the results de- 
rived from the 7-year WMAP data in combination with 
lower-redshift measurements of the expansion rate [13], 
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BBN gives smaller uncertainties in € by a factor of 4 
because there is nearly no correlation between the defor- 
mation parameter and the baryon density parameter as 
shown in Figure 2. The estimate of Q,h? is agreement 
with that from the CMB data [2] with errors. 

In summary, we have shown that the BBN puts strong 
constraint on the deformed parameter in the Lorentz- 
violating extension of the standard model, € = 0.036 + 
0.023. Since the BBN-predicted abundance of helium-4 
is very sensitive to the deformed parameter but less sen- 
sitive to the baryon energy density parameter, there is 
nearly no correlation between the two parameters. Our 
results indicate no significant preference for departure 
from Lorentz symmetry in the neutrino sector in the 
early Universe. Compared to previous constraints on 
the Lorentz-violating coefficient, current BBN observa- 
tion yields a weaker constraint. As listed in Table XIII 
of [12], the coefficient is constrained down to 107° from 
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